AN INTRINSIC APPROACH IN THE CURVED 
-BODY PROBLEM: THE NEGATIVE CURVATURE 

CASE 



Florin Diacu 
Pacific Institute for the Mathematical Sciences 

and 

Departament of Mathematics and Statistics 
University of Victoria 
P.O. Box 3060 STN CSC 
Victoria, BC, Canada, VSW 3R4 
diacu@uvic . ca 

Ernesto Perez-Chavela 
Departamento de Matematicas 
UAM-Iztapalapa 
Mexico, D.F. MEXICO 
epcOxanum . uam . mx 

J. Guadalupe Reyes Victoria 
Departamento de Matematicas 
UAM-Iztapalapa 
Mexico, D.F. MEXICO 
revg@xanum . uam . mx 

September 14, 2011 

Abstract. We consider the motion of n point particles of posi- 
tive masses that interact gravitationally on the 2-dimcnsional hy- 
perbolic sphere, which has negative constant Gaussian curvature. 
Using the stcrcographic projection, we derive the equations of mo- 
tion of this curved n-body problem in the Poincare disk, where we 
study the elliptic relative equilibria. Then we obtain the equations 
of motion in the Poincare upper half plane in order to analyze the 
hyperbolic and parabolic relative equilibria. Using techniques of 
Riemannian geometry, we characterize each of the above classes 
of periodic orbits. For n = 2 and n = 3 we recover some previ- 
ously known results and find new qualitative results about relative 
equilibria that were not apparent in an extrinsic setting. 
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1. Introduction 

We consider the negative curvature case of the curved n-body prob- 
lem, i.e. the motion of n point masses in spaces of constant Gaussian 
curvature under the influence of a gravitational potential that naturally 
extends Newton's law. This article is the completion of an earlier study 
done for positive curvature by Ernesto Perez-Chavela and J. Guadalupe 
Reyes Victoria, [26J. The novelty of the approach used in these twin 
papers is the introduction of intrinsic coordinates, which until recently 
seemed to be out of reach due to the elaborated computations they 
involve. But after Florin Diacu, Ernesto Perez-Chavela, and Manuele 
Santoprete obtained the general equations of motion for any number n 
of bodies in terms of extrinsic coordinates, j8] , a glimpse of hope ap- 
peared for an intrinsic approach, which has been now achieved for the 
2-dimensional case in [26] and the present paper. An intrinsic study of 
the 3-dimensional case is yet to be done. 

Whereas [26] had set a basic strategy, which we followed here too, in 
the present study we met with more technical difficulties. In [26], the 
idea to use stereographic projection and analyze the relative equilibria 
in the spherical plane had been successful, but the similar approach 
proved not strong enough for negative curvature. Apart from deriving 
the equations of motion in the Poincare disk to study elliptic relative 
equilibria, we also had to use the Poincare upper-half-plane model of 
hyperbolic geometry in order to overcome the hurdles encountered in 
the disk model when analyzing hyperbolic and parabolic relative equi- 
libria. This outcome seems to confirm an old saying, which claims that 
"in celestial mechanics, there is no ideal system of coordinates." The 
results are, of course, model independent, but it appears that each of 
them is easier to prove and understand with the help of a particular 
model. 

1.1. Motivation and history. The motivation behind the study of 
the curved n-body problem runs deep. In the early 1820s, Carl Frie- 
drich Gauss allegedly tried to determine the nature of the physical 
space within the framework of classical mechanics, i.e. under the as- 
sumption that space and time exist apriori and are independent of 
each other, [21], [J3]. He measured the angles of a triangle formed by 
three mountain peaks, hoping to learn whether space is hyperbolic or 
elliptic. But the results of his measurements did not deviate form the 
Euclidean space beyond the unavoidable measurement errors, so his 
experiments were inconclusive. Since we cannot reach distant stars to 
measure the angles of large triangles, Gauss's method is of no practical 
use for astronomic distances either. 
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But celestial mechanics can help us find a new approach towards 
establishing the geometric nature of the physical space. If we extend 
Newton's n-body problem beyond the Euclidian case and also prove 
the existence of solutions that are specific to each of the negative, zero, 
and positive constant Gaussian curvature spaces, then the problem 
of understanding the geometry of the universe reduces to finding, in 
nature, some of the orbits proved mathematically to exist. 

Therefore obtaining the natural extension of the Newtonian n-body 
problem to spaces of non-zero constant Gaussian curvature, and study- 
ing the system of differential equations thus derived, appears to be a 
worthy endeavour towards comprehending the geometry of the gravita- 
tional universe. Additionally, an investigation of this system when the 
curvature tends to zero may help us better understand the dynamics 
of the classical case, viewed as a particular problem within a general 
mathematical framework. 

The attempts to find a suitable extension of Newton's gravitational 
law to spaces of constant curvature started in the 1830s with the work 
of Janos Bolyai, [2], and Nikolai Lobachevsky, [22], who considered a 
2-body problem in hyperbolic space. Unfortunately, their purely geo- 
metric approach led to no immediate results. An important step for- 
ward was made in 1870 by Ernest Schering, who expressed the ideas of 
Bolyai and Lobachevsky in analytic terms, [21]. This achievement led 
to the formulation of the equations of motion of the 2-body problem, 
both for negative and positive curvature. Since then, many researchers 
brought contributions to the curved 2-body problem, including Wil- 
helm Killing, [TS], [TB], [T7], and Heinrich Liebmann, [Tj5], [20], [21] • 
The latter showed that the orbits of the curved Kepler problem (i.e. 
the motion of a body around a fixed center) are conies in the hyperbolic 
plane and proved an analogue of Bertrand's theorem, [I], which states 
that for the Kepler problem all bounded orbits are closed. But the 
most convincing argument that the extension of Newton's law due to 
Bolyai and Lobachevsky is natural rests with the fact that the poten- 
tial of the Kepler problem is a harmonic function in the 3-dimensional 
space, i.e. a solution of the Laplace- Beltrami equation, in analogy with 
the classical potential, which is a solution of Laplace's equation, i.e. it 
is a harmonic function in the classical sense, [IB] . 

Recently, Florin Diacu, Ernesto Perez-Chavela, and Manuele Santo- 
prete proposed a new setting for the problem, which allowed an easy 
derivation of the equations of motion for any n > 2 in terms of extrin- 
sic coordinates, [8]. The combination of two main ideas helped them 
achieve this goal: the use of Weierstrass's hyperboloid model of hyper- 
bolic geometry (also called the hyperbolic sphere) and the application 
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of the variational approach of constrained Lagrangian dynamics, [12] . 
They also succeeded to solve Saari's conjecture in the collinear case 
for both positive and negative curvature, having settled the Euclidean 
case earlier, [9], [7]. Remarkable is also their discovery of the fact that, 
in the curved 3-body problem, all Lagrangian orbits (i.e. rotating equi- 
lateral triangles) must have equal masses. Since Lagrangian orbits of 
non-equal masses exist in our solar system, such as the rotating equi- 
lateral triangles formed by the Sun, Jupiter, and the Trojan asteroids, 
it means that, at distances of the order of 10 1 AU, space is Euclidean. 

Since then, several papers analyzed the equations of motion. Florin 
Diacu studied the singularities of the equations and of the solutions, 
[5], Florin Diacu and Ernesto Perez-Chavela provided a complete clas- 
sification of the homographic orbits in the curved 3-body problem, 
[6], Florin Diacu considered the polygonal homographic orbits for any 
finite number of bodies, [I], and Pieter Tibboel solved an open prob- 
lem stated in [I]. Along the same lines, Regina Martinez and Carles 
Simo studied the stability of the Lagrangian relative equilibria and ho- 
mographic orbits for the unit sphere, [23| . All these papers treat the 
2-dimensional case. The only study, so far, of the 3-dimensional curved 
n-body problem is [3] , in which Florin Diacu analyzed relative equilib- 
ria, a paper in which more details on the history of the problem, as 
well as an extensive bibliography, can be found. 

1.2. Our results. As mentioned earlier, this paper is a natural contin- 
uation and completion of , which was a study of the curved n-body 
problem with the help of intrinsic coordinates in 2-dimensional spaces 
of positive curvature. Here we consider the negative curvature case. 

In section 2, we start with the extrinsic description of the motion of 
n point particles of positive masses on the hyperbolic sphere h 2 R , given 
by the upper sheet of the hyperboloid of 2 sheets, 

(1) x 2 + y 2 - z 2 = -R 2 , z>0, 

embedded in the Minkowski space Mf. Using stereographic projection 
through the north pole of the hyperbolic sphere ~h 2 R , we move the prob- 
lem to the Poincare disk D^, where we obtain the equations of motion 
in complex coordinates and show that the original equations of motion 
are equivalent to the equations obtained in D^. Then we provide the 
first integrals of the system in J} 2 R . 

In section |3j we consider a suitable Killing vector field in 3 R and 
its associated one-dimensional additive subgroup of the Lie algebra 
su(l, 1). When we project this Lie algebra via the exponential map 
onto the Lie group £77(1, 1), we obtain a one-dimensional subgroup of 
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proper isometries which generate the elliptic relative equilibria of the 
curved n-body problem in 3 2 R . By analyzing the corresponding one- 
dimensional subgroup of Mobius transformations, we obtain the equa- 
tions that characterize all the elliptic relative equilibria. Then we study 
these elliptic relative equilibria for n = 2 and n = 3. In the former case 
we give a complete description of the elliptic relative equilibria. In the 
latter case we study both the Eulerian and the Lagrangian orbits and 
describe their qualitative behavior, thus recovering in 3 2 R some results 
proved in [5] for L^. Additionally, we show that other elliptic rela- 
tive equilibria don't exists in the curved 3-body problem in the case of 
negative curvature. 

Although we also obtain in section 3 the one-dimensional subgroups 
of proper isometries that generate the hyperbolic and parabolic rela- 
tive equilibria, their complicated expressions make these orbits hard to 
understand. Therefore, in section HI we move the problem from the 
Poincare disk to the Poincare upper half plane M 2 R , in which we 
obtain the intrinsic equations of motion. Again, we consider the suit- 
able Killing vector fields in M 2 R and their corresponding one-parameter 
additive subgroups in its Lie Algebra s[(2, M), which we project via the 
exponential map onto the Lie group of proper isometries SL(2, M). The 
subgroups corresponding to elliptic relative equilibria become compli- 
cated in this context, but we already studied them in the Poincare disk. 
Fortunately, the subgroups corresponding to hyperbolic and parabolic 
relative equilibria are simple enough to allow a complete analysis. So 
using the associated one-parameter subgroups of Mobius transforma- 
tions, we obtain the conditions for the existence of hyperbolic and 
parabolic relative equilibria. Then we proceed as in section 3 with an 
analysis of the cases n = 2 and n = 3 for the former orbits, whose 
qualitative behavior we describe if they are of Eulerian type. Finally 
we confirm, within our more general context, a result obtained in [8], 
which proves that parabolic relative equilibria do not exist. 

2. Equations of motion 

Consider a connected and simply connected 2-dimensional surface of 
constant negative Gaussian curvature n = —1/R 2 . It is well known 
(see, e.g., [TU]) that this surface can be locally represented by ([T]), 
the upper sheet of the 2-dimensional hyperbolic sphere of radius iR, 
denoted by L^, embedded in the Minkowski space Rf , which is endowed 
with the Lorentz inner product defined by 

(2) Qi Q 2 := ?i?2 + Ditfe - 3i32 

for any points Qi = t)i,3i) and Q 2 = (? 2 , ^2,32) of K?. 



6 F. Diacu, E. Perez-Chavela, and J.G. Reyes Victoria 

Let II be the stereographic projection of through the north pole, 
(0,0, —R), located at the vertex of the lower sheet of the hyperbolic 
sphere. The image of II is the planar disk of radius R, denoted by 3 2 R , 
having the center at the origin of coordinate system, so we can write 
this function as 

(3) U:h 2 R ^B% II(Q) = g. 

It is easy to see that n(L^) is the entire open disk 3 R , and that if we 
take Q = (p, tj, 3), then q = (u, v), where 

(4) u = -^- and Rt) 



R + d R + d 

The inverse of the stereographic projection is given by the equations 



2R 2 u 2R 2 v R(R 2 + u 



J + v 2 

R 2 — u 2 — v 2 ' R 2 — u 2 — v 2 ' R 2 — u 2 — v 2 

The metric (distance) of the sphere h 2 R is transformed into the metric 

4R 4 

(6) ds2= (Ri-u 2 -v 2 ) 2 (du2 + dv2) - 

In the Poincare disk B> R , the geodesies are the diameters of the circle 
of radius R and the arcs of the circles orthogonal to it. In terms of the 
above metric, all geodesies have infinite length, so the circle of radius 
R represents the points at infinity. 

Since the metric is conformal, with factor of conformity 

AR A 

XM = (R 2 -u 2 -v 2 r 

the Christoffel symbols associated to the metric are given by 

1 d\ 1u 



~2X( Ul v)du (i? 2 — u 2 — v 2 ) 2 ' 
1 d\ 2v 



(7) - V\ 2 = r 2 n = r\ 

(8) T 2 22 = -T l n = Y\ 2 - ^-^^ - {R2 _ u2 _ v2r 

so the geodesies can also be obtained by solving the system of second 
order differential equations (see [TT] for more details): 

u + T^u 2 + 2T\ 2 uv + T\ 2 v 2 = 0, 

v + T 2 u ii 2 + 2T 2 12 uv + T 2 22 v 2 = 0, 
which is equivalent to 

,,> (uii 2 + 2vuv - uv 2 



R 2_l-2_ V 2 (vv 2 + 2uuv - vu 2 ) 
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From now on, we will think of 3 2 R as being endowed with the above 
metric. The geodesic distance between two points q k ,qj £ ^> R is 

Qk © Qj 



d k j = d(q k ,qj) = d(Q k ,Qj) = -Rcosh 1 ^~ 



R 2 



In analogy with [26] , where we considered the cotangent potential for 
k > 0, we will use here the hyperbolic cotangent potential for n < 0. 
With the help of the stereographic projection, we see that 

Qk Qj 



4R*q k - qj -R 2 (\\q k \ 


\ 2 + R 2 )(\\ qj \\ 2 + R 2 ) 


(R 2 -\ 


\qk\ 


\ 2 )(R 2 -H\\ 2 ) 



where • is the standard inner product in the plane M. 2 . Then 
(10) coth R / 4 ^ - ^ R2( lk-q j -{\\q k \\ 2 + R 2 ){{\\q j \\ 2 + R 2 

where 



R W 



w = ^/W x - w 2 , 

Wi = [AR 2 q k ■ qj - (\\q k \\ 2 + R 2 )((\ \ Qj \ \ 2 + R 
W 2 = (R 2 -\\q k \\ 2 mR 2 -\\ qj \\ 2 ) 2 . 
We now introduce complex variables in 3 2 R , 

z = u + iv, z = u — iv, 

in order to simplify the computations. The inverse of this transforma- 
tion is given by 

z+z z-z 
11 u = , V = . 

Then in terms of z and z, the formulas, (J5J), for the inverse of the 
stereographic projection are given by 

R 2 (z + z) %R 2 {z-z) R(\z\ 2 + R 2 ) 



(12) 



R 2 -\z\ 2 ' ' R 2 -\z\ 2 ' R 2 -\z\ 



where |.| denotes the absolute value of a complex number. 
The distance © and equation (fTDl take the form 

AR A 

(13) ds 2 = —z — pjrx dzdz 

{R z — \z\ z ) z 

and 

(14) coth fl 



2) 



dkA = 2(z k Zj + ZjZ k )R 2 - (\z k \ 2 + R 2 )(\zj\ 2 + R 

R J \f®2,(k,3){z,Z) 

where 

(15) Q2,W)(z,z) = [2(z k Zj + z jZk )R 2 - {\z k \ 2 + i? 2 )(|^| 2 + R 2 )} 2 
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-{R 2 -\ Zk \ 2 f{E?-\ Zj \ 2 f 

and dkj = d(z k ,Zj) denotes the geodesic distance in the metric (jT3l) 
between the points z k and Zj in 3 R . From now on we will think of 
B> R as the hyperbolic Poincare disk endowed with this new form of the 
metric. 

Notice that, in these complex coordinates, system which yields 
the geodesies, takes the form 

/ \ 2zi 2 

(16 » s + w^w = °- 

2.1. The intrinsic approach. Let z = (zi,Z2, ••• ,z n ) E (B#) n be 
the configuration of n point particles with masses mi,rri2, ■ ■ ■ , m n > 
in B> R . We assume that the particles are moving under the action of 
the Lagrangian 

(17) L R (z, z, z, z) = T R (z, z, z, z) + U R (z, z), 
where 

(18) T R (z, z, z, z) = - ^ m k \(z k ,z k )\z k \ 2 



2 

Kk<n 



is the kinetic energy and 
(19) U R (z,z) = 



l - J2 rn k m jCoih R f-^J 

l<fe<j<n 



n 

4 £ 



m k mj 



2(z k z i + ^ fc )i? 2 - ( 2 fc M + # 2 )(N 2 + R 2 



is the force function of the particle system {—Ur is the potential) de- 
fined in the set 3 2 R \ A, where A is the set of zeroes of 02,(fcj)( 
given by (IIS]) , and 

4R 4 

(20) A(z fc ,%) 



z, z 



(R 2 -\z k \ 2 ) 2 

is the conformal function of the Riemannian metric. 

The following result gives the equations of motion for the problem. 
It's proof is perfectly similar with the homologue result for curvature 
k > 0, obtained in 126 1, so we omit it here. 



Lemma 1. Let 

1 n 

L(z, z) = - ^2 m kK z k, z k )\z k \ 2 + U R (z, z) 

k=l 
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be the Lagrangian defined in ( TT71) for the given problem. Then the solu- 
tions of the corresponding Euler- Lagrange equations satisfy the system 
of second order differential equations 



2m k z k z 2 k 2 <9£7 



.it- 



(21) m k z k = -— + -J—, k = 1, . . . , 

ti — \z k \ A{z k ,z k ) oz k 

where 

, , 9Ur _ >A 2m k m j RP 2t(kJ) (z,z) 

[ j dz k -2-, (e 2 , {kJ) (z,z))^ ' 

P2,(k,j)(z, z) = {R 2 - \z k \ 2 ){R 2 - \z d \ 2 ) 2 {z d - - z k ){R 2 - z k z 

©2,(fcj)CM) = + %^)^ 2 - d^l 2 + R2 )(\ Z j\ 2 + R2 )T 

-(R 2 - \z k \ 2 ) 2 (R 2 - \ Z] \ 2 ) 2 , k, j e {1, ... , n}, k ^ j. 

Notice that the first term on the right hand side of ([21]) depends 
on the kinetic energy alone, whereas the second term depends on the 
potential. Therefore, from Lemma [I], we can also draw the following 
conclusion. 

Corollary 1. If in 3 2 R the potential is constant in the entire space, 
then the particles move freely along geodesies. 

Proof. Since the potential is constant, equations (l2Tj) take the form 

(23) m k z k + — j — - = 0, k = 1, . . . , n, 

R 1 - \zk\ 

and, if we simplify by m k > 0, we obtain that the coordinates of each 
body satisfy equation ( fToT) . so each body moves along a geodesic. □ 



2.2. Equivalence of the models. In [8], where V = (<9 r , d^, —d s ), the 
equations of motion for the n-body problem in the hyperbolic space 
are 

(24) mjQj = V Qj V K (Q) - rrijKiQj Q 3 )Q 3 , j = 1, . . . , n, 

where Q = (Qi, • • • , Q n ) is the configuration of the system, V K is the 
force function, 

(25) V K (Q)= £ rmmWKQtOQ, 



, j<n [(nQi Qj) 2 - (kQi Q 1 ){kQ 3 Qj)] 1 ' 2 ' 

Qj OQj = k _1 = R 2 , j = 1, . . . , n, are the constraints, which maintain 
the particles on the hyperbolic sphere, and denotes the Lorentz inner 
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product of the Minkowski space, denned in (jSJ). The (^-gradient, / 
1, . . . , n, of the force function is 



(26) v Q MQ)= E [9j " ( ; g 1 / ? j)g ' 1 

j=ij# [(kQi Qjf - l] 

As shown in [8J, system f !24p has 4 integrals of motion. One is the 
integral of energy, 

1 - 

(27) - J2 m AQj © Qj){*Qi © Qi) - v«(Q) = h, 

3=1 

where /i is the energy constant, and the other 3 are the integrals of the 
total angular momentum, 

n 

(28) ^mjfaij - iji),) = c u 

i=i 



(29) J2 m j(hh ~hh) = c ^ 

i=i 

(30) E m i(^?i-fiOi) =c 3 , 

i=i 

where C\,C2, C3 are the constants of the total angular momentum. 

Now we can state the equivalence between the equations of motion 
of the curved n-body problem defined on the hyperbolic sphere L^, 
of Gaussian curvature k = — 1/R 2 < 0, and the equations we have 
obtained in the Poincare disk D^. The proof of this result is perfectly 
similar with the proof of Theorem 2.3 in which concerns the case 
k > 0, so we omit it here. 

Theorem 1. The equations of motion of the curved n-body problem 
on the Poincare disk 3 2 R! given by system (1211) . and the correspond- 
ing equations on the hyperbolic sphere h 2 R , given by system f)24p . are 
equivalent. 

2.3. Integrals of motion. Since Theorem [1] proves the equivalence 
between the equations of motion of the curved n-body problem in 
and B) R , system (121!) inherits the properties of system (!2~4l) . We are in- 
terested here in how the first integrals of motion translate from system 
(|2"4"|) to system (|2T|) . Using equations (fT2"|) . which give the inverse of 
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the stereographic projection, we see that the integrals of motion given 
by (EZD, d2HD, flUD, and ([HD become, respectively, 



(31) ^^2m k X(z k ,z k )\z k \ 2 - U R (z,z) = h, 

k=l 

where h is the same energy constant as in (127|) . and 



(32) LpT^)2^ ~ Zk > + Zk ^ k ~ ZkZ M = Cl ' 



fc= 



(33) ^[i? 2 (i fc + 4) -4^] = c 2 , 
fc=1 I I k\ ) 

(34) ^ /■ B 2_u,l2^ ^ fcifc " Z * J * ) = Cs ' 



A straightforward computation confirms that the left hand sides of 
equations (1321) . fl33|) . and ( 134")) are real functions. Moreover, 01,02,03 G 
IR are the same constants as in ([28]), fl29l). and (1301). 



3. Relative equilibria in D| 

In this section we start to analyze the dynamics of the particles that 
interact in Our goal is to give some general characterization of the 
relative equilibria, which we define below with the help of geometric 
tools. For this purpose, we will need to understand first what are the 
singularities of the equations of motion, such that to avoid the singular 
configurations when dealing with regular equilibria, and to introduce 
the Principal Axis Theorem, which will guide us towards finding a 
proper definition of the relative equilibria. 

3.1. The singularities of the equations of motion in D^. The 

equations of motion (1211) have singularities in ©|, when at least one 
denominator vanishes, i.e. if there exist indices k, j, with 1 < k < j < 
n, such that 

< S>2,(k,j)(z,z) = 0. 

A straightforward computation shows that this statement is equivalent 
to saying that there are k, j, with 1 < k < j < n, such that 

(R 2 - z K Zj)(R 2 - ZjZ K )(z K - Zj){zj - z K ) = 0. 
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Since \zi\, \zi\ < R, i = 1, . . . , n, it follows that this equation is satisfied 
only for Zj = Zk, i.e. when a collision takes place. Therefore the set of 
singularities of the equations of motion is 

A = (J Afcj, where A kj = {(z u . . . , z n ) \ z k = Zj}. 

l<k<j<n 



3.2. Principal axis theorem in L|,. Let 5*0(1, 2) denote the Lorentz 
group, defined as the Lie group of all the isometric rotations of deter- 
minant 1 in that keep the hyperbolic sphere h R invariant. An im- 
portant result related to this subgroup of orthogonal transformations 
is the Principal Axis Theorem, [25], which states that there are 3 1- 
parameter subgroups of S0{1,2) whose elements can be represented, 
in some suitable basis of R?, as 



(1) matrices of the form 



(35) 



A = P 



cos 9 — sin 9 
sin 9 cos 9 
1 



P- 



called elliptic rotations in h 2 R around a timelike axis, the 3-axis 
in our case, 
(2) matrices of the form 



(36) 



B = P 



1 
cosh s sinh s 
sinh s cosh s 



called hyperbolic rotations in h 2 R , around a spacelike axis, the 
y-axis in our case, and 
(3) matrices of the form 



(37) 



C = P 



1 -t t 
t l-t 2 /2 t 2 /2 
t -t 2 /2 l + t 2 /2 



P- 



called parabolic rotations in L|,, around a lightlike axis, which 
is given here by j: = 0, t) = 3, where P G SO(l, 2) in all cases. 

Remark 1. This result implies that any isometry in h R can be written 
as a composition of an elliptic rotation around the 3-axis, a hyperbolic 
rotation around the jr-axis, and a parabolic rotation around the axis 
t) = 0,3 = r. 



An intrinsic approach in the curved n-body problem of negative curvature 



13 



3.3. Relative equilibria in D^. Let Iso(D^) be the group of isome- 
tries of Ofj, and assume that {G(t)} is a 1-parameter subgroup of 
Iso(D|j) that leaves 3 2 R \ A and A invariant. We can now give a 
general definition for relative equilibria in 3 2 R . 

Definition 1. A relative equilibrium of the negatively curved n-body 
problem is a solution z of equations (|2"T|) that is invariant relative to 
some subgroup {G{t)} oflso{D 2 R ), i.e. the function w, given by w{t) = 
G{t)z{t), obtained by the action of some element ofG, is also a solution 
of system (T2TT) • 

To understand the implications of this definition and be able to rep- 
resent the relative equilibria in 1D) R in as precise terms as we did for h 2 R 
in [8], we need to take first a look at the topological group structure 
of isometric rotations of 3 2 R (for more details, see, e.g., [H]). For this 
purpose, consider the matrix 



be the special orthochronous unitary group. Then some algebraic com- 
putations show that any matrix A G SU(1, 1) has the form 



with a, b G C satisfying \a\ 2 — \b\ 2 = 1. This last condition implies that 
the group SU(1, 1) is diffeomorphic with the real 3-dimensional unit 
hyperbolic sphere embedded in C 2 . The term orthochronous means 
that the transformations A do not change the direction of the time t 
in the standard interpretation of the Minkowski space. In our case, 
however, this is just another space coordinate. 

The group of proper orthochronous isometries of D^, i.e. the group 
of transformations that also maintain the geometric orientation, is the 
quotient group 



where I is the unit 2x2 matrix. To every class A G SU(1, l)/{±/}, 
we can associate a Mobius transformation, 




and let 



SU(1, 1) = {A G GL(2, C) | A T IA = /, det A = 1} 




SU(l,l)/{±/}, 



f A :B 2 R ^B 2 R , f A (z) 



az + b 



bz + a' 



for which it is easy to see that /-a(^) = /a(^)- 
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If M(2, C) is the set of 2 x 2 matrices with complex elements, the 
Lie algebra of SU(1, 1) is the 3-dimensional real linear space 

su(l,l) = {X e M(2,C) | IX T = -XI, trace X = 0} 

spanned by the Killing vector fields in H) 2 R associated to the Pauli ma- 
trices, 

1 / 1 \ 1 / i \ 1/0 i 



91 - 2 v i o ; ' 92 - 2 o -i ; ' ^ 3 - 2 v 

which form a basis of su(l, 1). 

Consider further the exponential map of matrices, 

exp: su(l,l) -> SU(1,1), 

applied to the one-parameter additive subgroups {tg\}, {tg2}, and 
{i<?3}, which are straight lines form the geometric point of view. This 
operation leads us to the following one-parameter subgroups of SU(1, 1) 
in D|: 

(1) the subgroup 

C (t) ~ nm(ta ) - ( cosh ( t / 2 ) sinh(t/2) \ 
G 1 (t) - expitgi) ~ { sinh(t/2 ) cosh(t/2) ) ' 

which defines the one-parameter family of Mobius transforma- 
tions 

cosh(t/2)2: + sinh(t/2) 
m Ml ' j "sinh(t/2)^ + cosh(t/2)' 

(2) the subgroup 

e «/2 o 



G 2 (t) = exp(tg 2 ) 



e^ 2 



which defines the one-parameter family of Mobius transforma- 
tions 

(39) f G2 (z,t) = e u z, 

(3) the subgroup 

C (t) - nroffo 1 - ( cosh (*/ 2 ) *sinh(i/2) \ 
G 3 (t) - exp{tg 3 ) - ^ _ zsinh(t/2 ) cosh(t/2) ) ' 

which defines the one-parameter family of Mobius transforma- 
tions 

cosh(t/2)2 -Msinh(t/2) 



(40) f G3 (z,t) 



-i sinh(t/2)z + cosh(t/2)' 
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3.4. Elliptic relative equilibria. When projected into the Poincare 
disk, a compositions of the elliptic, hyperbolic, and parabolic rotations 
of h 2 R corresponds to some composition of the Mobius transformations 
( )38|) . (|39|) . and (|40|) . We therefore need to analyze each of these one- 
parameter subgroups of transformations in order to find their associ- 
ated relative equilibria. We will first achieve this goal for the second 
group, G2, given by the Mobius transformations /g 2 (- z ) = eltz i which 
correspond to the elliptic rotations in L^. The other 2 groups will 
be treated in section 4, since their study becomes very tedious in B) 2 R 
because of complicated computations. 

Notice that the one-parameter rotation subgroup of 50(1,2) intro- 
duced by the matrix 



(41) A(t) 



i.e. the rotation (|35|) around the 3 axis of that leaves invariant any 
hyperbolic sphere h R centered at the origin, is the isometric flow for 
the basic Killing vector field 

(42) M?>«u) = H),?,o). 

This observation leads us to the following result. 

Proposition 1. LetH: SU(1, l)/{±/} — > S0(l,2) be an isomorphism 
between the groups of proper orthochronous isometries of the Poincare 
disk H) R and the Lorentz group of the hyperbolic sphere h R . Then 

H(G 2 (t))=A(t). 
Proof. The stereographic projection, 

shows that since the rotation tangent vector at (y, rj, 3) in h 2 R is (— t), y, 0), 

after the projection we have 

(43) 

/ _R_ R " , \ f -*)\ ( -M. 

Dn[(L z )( P , n ,3)] r = y R W ? = R A 

V U R +i (R+i) 2 J \ J \ R+i 

where DU is the Jacobian matrix and the upper T denotes the trans- 
posed of the vector. In complex notation this relationship corresponds 
to 

— v + iu = i(u + iv) = iz, 
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which leads us to the differential equation 

(44) z = iz. 

Its flow is given by ft(z) = e lt z, associated to the one-parameter 
subgroup of Mobius transformations fa 2 - This remark completes the 
proof. □ 

In order to obtain from Proposition [1] some information regarding 
the relative equilibria of type G2 in D^, we consider functions of the 
form 

(45) w k (t) = e u z k (t), 

where z = (zi, . . . , z n ) is a solution of equation (1211) . and look for 
conditions that the function w = (wi, . . . ,w n ) is also a solution of 
system (T2"Tj) . Straightforward computations show that 



(46) 



= [iz k + z k )e lt 
= (z k + 2iz k - z k )e lt 
#- = e ft , k = l n. 



Using these facts together with the conditions that w is a solution of 
equation fl2Tj) . 



m k w k 



2m k w k w\ (R 2 - \w k 



2\2 



R 2 



w k 



AR A 



du A k 



dwi 



we obtain in terms of z that 

m k Vz k + 2iz k 



2m k e- u z k {iz k + z k ) 2 e 2U {R?_ 



Zk 



2\2 



R 2 



Zk 



2m k z k (iz k + z k Ye 



2„it 



(R 2 



AR 4 

\zk 



9Ur dz k 
dz k dw k 



2\2 



dU. 



R 2 



R e u . 



Zk 



4i? 4 dz k 

Since z is a solution of ( }2"T|) . m k ^ 0, and e lt 7^ 0, the last relationship 
becomes 

2z k (2iz k z k — zl) 



(47) 



2iz k - z k 



R 2 



Zk 



which is equivalent to the equation 



(48) 



2% 



1 + 



2 


z k 


2 


R 2 - 




Zk 


2 



Zk 



1 + 



2 


Zk 


2 


R 2 - 




Zk 


2 



Zk- 



Equation ( )48|) holds if and only if 



1 + 



2 


Zk 


2 


R 2 - | 


Zk 


2 



or 2iz k = z k , k = l,...,n. 
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The first set of conditions is equivalent to |zfc(£)| = = R, which 
never holds. The second set of conditions, which provides information 
about how the velocities must behave for this kind of relative equilibria, 
holds for = r k , where r k > 0, k = 1, . . . , n. For this reason, the 

particles form a relative equilibrium associated to the Killing vector 
field (142|) if they are moving along Euclidean circles centered at the 
origin of the coordinate system in D|j. In terms of "L 2 R , the particles 
move on circles obtained by slicing the hyperbolic sphere with planes 
orthogonal to the vertical axis, 3, of M.\. 
We can now prove the following result. 

Theorem 2. Consider n point particles with masses m 1; . . . , m n > 0, 
n > 2, moving in 3 2 R . A necessary and sufficient condition for the 
function z = (z±, . . . , z n ) to be a solution of system f[2Tj) that is a relative 
equilibrium associated to the Killing vector field Lz defined by equation 
( 142|) is that for all k — 1, . . . , n, the following equations are satisfied 

(AO) R 3 (R 2 + r 2 k )z k _ A m 3 [r) - R*)\ Zj - z k ){R? - z k z 3 ) 
where r k = \z k \ and 

@2,(k,j)(z, z) 

= [2{z kZj + zfz k )R* - (rl + i? 2 )(r? + i? 2 )] 2 - (i? 2 - r 2 ) 2 (i? 2 - r 2 ) 2 , 
k,j e {1, . . . ,n}, k^ j. 

Proof. From the equations 2iz k = z k , k — 1, . . . , n, we concluded that 
a necessary condition for the existence of a relative equilibrium of the 
aforementioned type is that the particles move along ordinary circles 
centered at the origin of the coordinate system in D^. Differentiating 
these conditions and using them again we obtain that 

(50) - A'z k = z k , k = l,...,n. 

Comparing these equalities with the equations of motion (12T1) . we con- 
clude that the coordinates of a relative equilibrium must satisfy the n 
algebraic equations 



(51) m k z k = 



2m k \z k 


\ 2 z k 2(R 2 - 


Zk 


2 ) 2 dU R 


R 2 - 






2 i? 4 


dz k 



, k = 1, . . . , n. 



Substituting r k = \z k \, k = l,...,n, into the above equations, we 
obtain the system of n equations ( |4"9]) . which characterize the relative 
equilibria given by the group G2- This remark completes the proof. □ 

Definition 2. We call elliptic relative equilibria the solutions of system 
( I2T|) that satisfy the conditions (149]) . 
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3.5. The case n — 2. We will next prove the existence of elliptic 
relative equilibria for 2 particles in H) 2 R , both in the case when they 
move on the same suitable circle and in the case when they move on 
different suitable circles. To achieve this goal, notice first that for n = 2 
and rrii, m2 > 0, equations (1491) take the form 

R 3 (R 2 + rj)z 1 m 2 (r 2 -R 2 ) 2 (z 2 -z 1 )(R 2 -z 1 z 2 ) 

m 2 -r\Y = [Q2,(i, 2) (^)] 3/2 

R 3 (R 2 + r 2 )z 2 mi (r 2 - R 2 ) 2 ( Zl - z 2 )(R 2 - z 2 z x ) 



m 2 ~rlY [Q 2 ,(2,i)(^)] 3/2 

where 

Q2,(k,j)(z, z) = [2^ 1 ^+ 22 ^ 1 ) J R 2 -(r?+ J R 2 )(r 2 2 + J R 2 )] 2 -(r?- J R 2 ) 2 (r 2 2 - J R 2 ) 2 . 

Some algebraic manipulations lead us to the equation 

(R 2 + r 2 )(R 2 - r 2 ) 2 m 1 _ z 2 (z 2 - Zl )(R 2 - Zl z 2 ) 
(R 2 + r 2 2 )(R 2 - r 2 ) 2 m 2 Zl ( Zl - z 2 )(R 2 - z 2 z x ) ' 

and if we simplify the right hand side we obtain 

( , (R 2 + rp(R 2 - rffmi R 2 z 2 - z x r 2 

[ ' (R 2 + r 2 )(R 2 - r 2 ) 2 m 2 R 2 Zl -z 2 r 2 ' 

This equation shows that there are no elliptic relative equilibria for the 
2-body problem in 3 2 R when one particle is fixed at the origin of the 
coordinate system. Indeed, if z\ =0, then n = 0, and the denominator 
of the right hand side vanishes. If z 2 = 0, then r 2 = 0, and the above 
equation becomes 

(R 2 + rl)R 2 mi _ 
(R 2 — ri) 2 m 2 

which has, obviously, no solutions. 

Equation (I52j) holds for all time, in particular when the particle mi 
reaches the real line. At that time instant, we have Z \ = a G M, and 
let us denote z := z 2 and r := \z 2 \. Then, if we solve equation ( I52p for 
z, we have 



rn^R 2 + a 2 ){R 2 - r 2 ) 2 R 2 - m 2 r 2 [R 2 + r 2 )[R 2 - a 2 ) 2 



a, 



m^R 2 + a 2 )(R 2 - r 2 ) 2 a 2 - m 2 R 2 (R 2 + r 2 )(R 2 - a 2 ) 2 

therefore z is also a real number, so either z = r or z = —r. In other 
words, when mi reaches the real line, m 2 reaches it too. Consequently 
the above equation becomes 

mi (R 2 + a 2 )(R 2 - r 2 ) 2 (±r)[R 2 - a(±r)} 



(53) 



m 2 (R 2 + r 2 )(R 2 -a 2 ) 2 a[R 2 - a(±r)] 
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Since a, r < R, it follows that the left hand side is positive, so the right 
hand side must be positive too. Consequently z = — r, and thus the 
above equation takes the form 

mi (R 2 + a 2 )(R 2 - r 2 ) 2 _ r 

^ 5 ^ m 2 {R 2 + r 2 ){R 2 - a 2 ) 2 ~ a' 

which we can use to estimate the value of r that gives the position of 
m 2 as a function of mi,m 2 ,R, and a. For this purpose, we consider 
the real function /: (-R, R) — > R, 

(55) f(x) = m ia (R 2 + a 2 )(R 2 - x 2 ) 2 - m 2 x(R 2 + x 2 )(R 2 - a 2 ) 2 , 

whose zeroes give us the desired values of r and, therefore, the elliptic 
relative equilibria for the 2-body problem in O^. For this purpose, let 
us first prove the following result. 

Lemma 2. The function f defined in (155j) has no double roots. 
Proof. The first derivative of the function / is 

(56) f\x) = -A mi ax(R 2 + a 2 )(R 2 - x 2 ) - m 2 (R 2 + 3x 2 )(,R 2 - a 2 ) 2 . 

Since i^O and —R < x < R, the double zeroes of / must satisfy the 
equations 

-Ax f{x) = 
(R 2 - x 2 ) f'{x) = 0, 
which are equivalent to the system 

A mi ax(R 2 + a 2 )(R 2 - x 2 ) 2 - Am 2 (R 2 x 2 + x A )(R 2 - a 2 ) 2 = 
-4 mi ax(R 2 + a 2 )(R 2 - x 2 ) 2 - m 2 [R 2 + 3x 2 )(R 2 - a 2 ) 2 (R 2 - x 2 ) = 0. 

If we add the equations of the above system, we obtain the quartic 
equation 

x 4 + 6R 2 x 2 + R A = 0, 

which has only non-real roots, ± V / -3±2 v / 2 R, a fact that completes 
the proof. □ 

We can now state and prove the following result, which characterizes 
the elliptic relative equilibria of the curved 2-body problem, expressed 
in terms of the Poincare disk model, D^, of the hyperbolic plane. 

Theorem 3. Consider 2 point particles of masses mi, m 2 > moving 
in the Poincare disk3 2 R , whose center is the origin, 0, of the coordinate 
system. Then a function z = (zi,z 2 ) is an elliptic relative equilibrium 
of system (l2Tj) with n = 2, if and only if for every circle centered at 
of radius a, with < a < R, along which mi moves, there is a unique 
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circle centered at of radius r, which satisfies < r < R and ( 1541) . 
along which m 2 moves, such that, at every time instant, mi and m 2 
are on some diameter of O 2 ^ with between them. Moreover, 

(1) if m 2 > mi > and a are given, then r < a; 

(2) if mi = m 2 > and a are given, then r = a; 

(3) if mi > m 2 > and a are given, then r > a. 

Proof. If the solution z of system ( 12 ip with n = 2 is an elliptic relative 
equilibrium, then equations ( H9|) are satisfied, and they lead to equation 
(I54p . from which we can compute r for given mi,m 2 > and a, with 

< a < R. Then, by Lemma [21 there is a unique r as desired, so the 
particles move as described, and the implication follows. 

To prove the converse, assume that for given mi,m 2 > and a, 
with < a < R, there is a unique r, which satisfies < r < R and 
(1M|) . such that the bodies move as described. Then the motion must 
be given by the function z = (zi,z 2 ) with 

Zi(t) = a(cost + i sint), z 2 {t) = — r(cost + i smt), 

with the relationship between mi,m 2 , a, and r given by (1511) . A straight- 
forward computation shows that this function is a solution of system 
(I2T|) with n = 2 and satisfies equations (|49l) . 

To see how the relative values of mi and m 2 determine the relation- 
ship between r and a, we evaluate the function /, defined in (j551) . at 
x = 0, a, R, and obtain 

/(0) = am^i? 2 + a 2 )R 4 > 0, 

/(a) = a{mi - m 2 )(R 2 + a 2 )(R 2 - a 2 ) 2 , 
f(R) = -2R 3 m 2 (R 2 - a 2 ) 2 < 0. 

Since, by (|56|) . the derivative /' is negative in the entire interval (0, R), 
the function / is strictly decreasing. The conclusion follows then from 
the above relationships. This remark completes the proof. □ 

3.6. The case n = 3. In the case of 3 particles in the Poincare disk 
H) R , equations (149]) become 

R\R 2 + r 2 )zi 

1 ' 4(i? 2 -r 2 ) 4 

m 2 {R 2 -r 2 ) 2 {z 2 -Zi){R 2 -ZiZ 2 ) 

{[2(ziz 2 + z 2 Zi)R 2 - (r 2 + R 2 )(r 2 + R 2 )] 2 - (R 2 - r 2 ) 2 (R 2 - r 2 ) 2 } 3 / 2 

m 3 (R 2 -r 2 ) 2 (z 3 -Zi)(R 2 -ZiZ 3 ) 

{[2(ziz 3 + z 3 Zi)R 2 - (r 2 + R 2 )(r 2 + R 2 )} 2 - (R 2 - r 2 ) 2 (R 2 - r 2 ) 2 } 3 / 2 ' 
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R\R 2 + r 2 )z 2 
[M) 4(R 2 - rlY 

m 1 (R 2 -r 2 1 )*(z 1 -z 2 )(R*-z 2 z 1 ) 

{[2(* 2 *i + z x z 2 )R 2 - (r 2 + R 2 )(r 2 + R 2 )} 2 - (R 2 - r 2 ) 2 (R 2 - r 2 ) 2 } 3 / 2 
m 3 (R 2 -r 2 3 ) 2 (z 3 -z 2 )(R 2 -z 2 z 3 ) 



(59) 



{[2(z 2 z 3 + z 3 z 2 )R 2 - {r 2 + R 2 ){r 2 + R 2 )} 2 - {R 2 - r 2 ) 2 {R 2 - r 2 ) 2 } 3 / 2 ' 

R\R 2 + r 2 )z 3 



A(R 2 - r 2 3 Y 



{[2(^3^ + Zl z 3 )R 2 - (r| + R 2 )(r 2 + R 2 )] 2 - (R 2 - r 2 ) 2 (R 2 - r 2 ) 2 } 3 / 2 

m 2 {R 2 -r 2 ) 2 {z 2 -z 3 ){R 2 -z 3 z 2 ) 

'{[2{z 3 z 2 + z 2 z 3 )R 2 - (r| + R 2 )(r 2 + R 2 )} 2 - (R 2 - r 2 ) 2 (R 2 - r 2 ) 2 } 3 / 2 ' 



3.6.1. Eulerian Solutions. We will start the study of the case n = 3 
with the Eulerian elliptic relative equilibria in for which the bodies 
lie on a rotating geodesic. Of course, we can assume that this geodesic 
rotates around the origin of the coordinate system, so then it must be 
a rotating diameter of D^. 

Assume that the particle mi is located at the center of the disk, i.e. 
z\ = 0, and that m 2 reaches at some time instant the positive axis of 
the real line, i.e. z 2 = r 2 =: a > 0. For m 3 , we denote z := z 3 and take 
|z 3 | = r 3 . Then equations (1571) . (|58|) . and (1591) become, respectively, 

Q = m 2 (R 2 -a 2 ) 2 _ m 3 z(R 2 -r 2 ) 2 ^ 
a 2 r? 



(61) 



R 3 (R 2 + a 2 )a m x a 



(62) 



4{R 2 - a 2 Y [(a 2 + R 2 ) 2 - {R 2 - a 2 ) 2 ] 3 / 2 

m 3 (R 2 -r 2 ) 2 (z-a)(R 2 - az) 

{[2a{z + z)R 2 - {a 2 + R 2 )(r 2 + R 2 )} 2 - (R 2 - a 2 ) 2 {R 2 - r 2 ) 2 } 3 / 2 

R 3 (R 2 + rl)z _ m\z 



A(R 2 - r 2 Y [{r 2 + R 2 ) 2 - (R 2 - r 2 ) 2 } 3 / 2 
m 2 (R 2 - a 2 ) 2 (a - z)(R 2 - az) 



{[2a(z + z)R 2 - (rf + R 2 )(a 2 + R 2 )} 2 - (R 2 - r 2 ) 2 (R 2 - a 2 ) 2 } 3 / 2 ' 

The following result will show that, when one particle is fixed at the 
origin of the coordinate system at the center of D^, there is just one 
class of Eulerian elliptic relative equilibria, namely orbits for which the 
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distance from the fixed body to the 2 rotating bodies is the same, and 
consequently those masses must be equal. In terms of the hyperbolic 
sphere h 2 R , the configuration is an isosceles triangle that rotates around 
its vertical height. 

Theorem 4. Consider 3 point particles of masses m 1 ,m 2 ,m 3 > 
moving in the Poincare disk 3 2 R , whose center is the origin, 0, of the 
coordinate system. Take a function z = (^1,^2,^3) that describes the 
positions of the particles, with Z\(t) = for allt. Then z is an Eulerian 
elliptic relative equilibrium of system (I2ip with n = 3 if and only if 
m 2 and m 3 are at the opposite sides of the same uniformly rotating 
diameter of a circle of radius a in H R , centered at 0, with < a < R, 
and mi = m 2 . 

Proof. From equation (I60l . we can conclude that z must be a negative 
real number, so z = —r 3 =: — r. This implies that (160]) becomes 

(63) 



a 2 r 2 



Let us further consider the previous equation together with (IBTj) and 
( 162|) in which we substitute the previous values of zi, r 1; z 2 , r 2 , z 3 , r 3 . 
Then we obtain the new system 



(64) 



(65) 



m 2 (R 2 — a 2 ) 2 m 3 (R 2 — r 



2\2 



a 2 r 2 



R 3 (R 2 + a 2 )a m x 



4{R 2 - a 2 Y 2 3 / 2 8i? 3 a 2 
m 3 (R 2 - r 2 ) 2 (r + a)(R 2 + otr) 



(66) 



{\-AR 2 ra - {a 2 + R 2 )(r 2 + R 2 )} 2 - (R 2 - a 2 ) 2 (R 2 - r 2 ) 2 } 3 / 2 ' 
R 3 (R 2 + r 2 )r mi 



A{R 2 - r 2 ) 4 " 2 3 / 2 8R 3 r 2 
m 2 (R 2 - a 2 ) 2 (a + r)(R 2 + ar) 



{[-4R 2 ra - (r 2 + R 2 )(a 2 + R 2 )] 2 - (R 2 - r 2 ) 2 (R 2 - a 2 ) 2 } 3 / 2 ' 

If we multiply both sides of (l6"5j) by m 2 (R 2 — a 2 ) 2 , and both sides of 
( l66l) by —m 3 (R 2 — r 2 ) 2 , when we add the resulting equations we get 

m 2 (R 2 + a 2 )a m 3 (R 2 + r 2 )r 



^ 67) (R 2 -a 2 ) 2 (R 2 -r 2 ) 2 
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From equations ( 1631) and ( 1671) we obtain the linear system having the 
masses m 2 and m 3 as unknowns, 



(6 



m 2 r 2 (R 2 - a 2 ) 2 - m 3 a 2 (R 2 - r 2 ) 2 = 
m 2 (R 2 + a 2 )(R 2 - r 2 ) 2 a - m^R 2 + r 2 )(R 2 - a 2 ) 2 r = 0, 

which has nontrivial solutions if and only if the principal determinant 
vanishes. A straightforward computations shows that this condition is 
equivalent to 

(R 2 + r 2 )(R 2 - « 2 )V - (R 2 + a 2 )(R 2 - r 2 ) V = 0, 

so the principal determinant vanishes when 

(R 2 + r 2 y = (R 2 + a 2 )a 3 
[ ' {R 2 - r 2 ) 4 (R 2 - a 2 ) 4 ' 

To find the values of r that solve equation ( 1691) . we consider the function 

(R 2 + x 2 )x 3 

(70) 9 :[0,R)^R, g(x) = { {R2 + _ J )4 , 

which is strictly increasing in its domain. Therefore equation ( 1691) holds 
only for x = r = a. If we substitute these values in the first equation 
of system (1681) . it follows that m 2 = m 3 . □ 

In Theorem HJ we required that all masses are positive. Let us now 
consider the case when mi = 0. We then obtain the following result. 

Proposition 2. Consider 3 point particles of masses m\ = and 
m,2 = WI3 > moving in the Poincare disk 3 2 R , whose center is the 
origin, 0, of the coordinate system. Take a function z = (z\, z 2 , Z3) 
that describes the positions of the particles, and assume that Zi(0) = 
and the real parts of z 2 {0) and 23(0) are 0, i.e. mi is initially at 
the center and m 2 ,m 3 are initially on the horizontal diameter ofD 2 R . 
Then a necessary condition for the particles to form an elliptic relative 
equilibrium, is that m 2 and m 3 rotate on the same suitable circle, being 
at every time instant at the opposite sides of some diameter of that 
circle, and m x lies at the the center of the disk for all time. 

Proof. Without loss of generality, we can take 

z x =r\ = z\ =: c, z 2 = 7-2 = z 2 =: a, and z 3 = z 3 = —f3. 

Then some straightforward computations show that equations (160|) . 
(16T|) . and (1621) become, respectively, 

R 6 (R 2 + c 2 )c m 2 (R 2 -a 2 ) 2 m 3 (R 2 - (3 2 ) 2 

\' ) Try?. „9\4 _ 0/ Z „\9/ D5 Z „\2 



(R 2 -c 2 Y 2(a - c) 2 (R 2 - ac) 2 2(/3 + c) 2 (R 2 + (3c) 2 ' 
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R\R 2 + a 2 )a m 3 (R 2 - (3 2 ) 2 



(72) 



{R 2 -a 2 ) 4 2(P + a) 2 (R 2 + (3a) 



2\2 



R\R 2 + f3 2 )f3 m 2 (R 2 ~a 2 ) 



(73) {R 2 -P 2 Y 2{(3 + a) 2 {R 2 + (3a) 2 ' 

Since m 2 = m 3 , equations f[72"j) and (175]) lead us to the relationship 

f74 x (fi 2 + a 2 )a = (i? 2 + /3 2 )/3 

1 J (i? 2 -a 2 ) 2 (i? 2 -^ 2 ) 2 ' 

Consider the smooth real function 



h: [0,R) -> E, 



(i? 2 + x 



2> 



(i? 2 - X 2 ) 2 



It is easy to see that h is strictly increasing, which implies that relation 
( 1731) holds if and only if a = (3. Therefore m 2 and m 3 must be at 
opposite sides of the same circle. 

If we now use the fact that a = f3 in equation flTTj) . its right hand 
side vanishes, so c = 0, a remark that completes the proof. □ 



3.6.2. Lagrangian Solutions. We will next study the Lagrangian ellip- 
tic relative equilibria in B> 2 R . As we proved in [8], such orbits, formed 
by rotating equilateral triangles, i.e. r := \zi\ = \z 2 \ = \z 3 \, exist only 
when the 3 positive masses are equal, m := mi = m 2 = m 3 > 0. The 
converse is also true, and Florin Diacu gave a proof of this result for 
polygons with n > 3 sides in [4j in the more general case of homo- 
graphic orbits, which allow not only rotation but also expansion and 
contraction of the configuration, from which the result for relative equi- 
libria follows as a particular case. Moreover, Pieter Tibboel recently 
proved that irregular polygons cannot form homographic orbits in the 
2-dimensional hyperbolic space, as well as on the 2-dimensional sphere 
as long as the motion is not along a great circle of the sphere, [28] . 
when the orbit is necessarily a relative equilibrium, a case in which 
non-equilateral triangles exist for suitable non-equal masses, as proved 
in [1] for the curved 3-body problem. The proof we give here for the 
case n = 3 follows the idea used for Theorem 5.3 in [26] for positive 
curvature. 

Theorem 5. Assume that 3 point particles of equal masses move along 
a circle of radius r centered at the center of the the Poincare disk D 2 R . 
Then a necessary and sufficient condition for the existence of an elliptic 
relative equilibrium is that the particles form an equilateral triangle. 
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Proof. With the values Z\ =: r, Z2 ='■ re 102 and z 3 = Z2 =: re* 6 * 3 , some 
straightforward computations bring equations (J57j) . (JSED, and flBUj) . re- 
spectively, to 

^(f^ + r 2 ) 



(75) 



4m(i? 2 - r 2 ) 6 
{e w * - 1)(R 2 -r 2 e- ie *] 



(76) 



{[4R 2 r 2 cos 9 2 - (r 2 + i? 2 ) 2 ] 2 - (R 2 - r 2 ) 4 } 3 / 2 
(e^ 3 - 1)(R 2 -r 2 e-^ 3 ) 
~{[4i? 2 r 2 cos# 3 - (r 2 + i? 2 ) 2 ] 2 - (i? 2 - r 2 ) 4 } 3 / 2 ' 

R 3 (R 2 + r 2 ) 
Am(R 2 - r 2 ) 6 

- - {[ 4jR 2 r 2 cos ^ _ ( r 2 + #2)2)]2 _ (^2 _ r 2)4}3/2 

e - ie2 (e^ 3 - e i02 )(i? 2 - rV^-^) 
~ {[4i? 2 r 2 cos(# 2 - 9 3 ) - (r 2 + R 2 ) 2 )} 2 - (R 2 - r 2 ) 4 } 3 / 2 ' 

fl 3 (fl 2 + r 2 ) 
4m(i? 2 - r 2 ) 6 

e -^ 3 (l -e^ 3 )(i? 2 -rV e3 ) 

- _ { [ 4jR 2 r 2 cog Q 3 _( r 2 + R 2y}2- {B? - r 2)4}3/2 

e -ie 3 ^ e ie 2 _ e ifeWjj2 _ r 2 e i(e 3 -6 2 )^ 

~ \[AR 2 r 2 cos(9 3 - 9 2 ) - (r 2 + i? 2 ) 2 )] 2 - (R 2 - r 2 ) 4 ] 3 / 2 ' 
Adding equations (1731) and (1761) . subtracting equation (1771) from the 
sum, and separating the real and imaginary parts of the resulting equa- 
tion, we obtain 

l-cos# 2 l-cos(#3-# 2 ) 2(1- cos 9 3 ) 



(77) 



(7 



D\2 D23 D13 



(79) + Sin(g3 ~ ° 2) = 0, 

-D12 -D23 

where 

D 12 = 8 3 / 2 i? 3 r 3 [l - cos 9 2 ] 3/2 [R 4 + r 4 - 2i?V cos # 2 ] 3 / 2 , 

D 23 = 8 3/2 i? 3 r 3 [l - cos(# 3 - e 2 )] 3/2 [R A + r 4 - 2i? 2 r 2 cos(fl 3 - 2 )f\ 

D 13 = 8 3 / 2 i? 3 r 3 [l - cos 9 3 f 2 [R A + r 4 - 2R 2 r 2 cos £ 3 ] 3/2 - 

It is easy to check that 9 2 = y and 9 3 = satisfy equations (178!) and 
therefore the configuration corresponds to an equilateral triangle. 
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Using standard trigonometry, equation ( 1791) becomes 



sin 2 (^; 
sin 2 (§) 



(r 2 - R 2 ) 2 + R 2 r 2 sin 2 (^; 
(r 2 - R 2 ) 2 + R 2 r 2 sin 2 (f) 



l_ sin 2(^; 
l-sin 2 (f) 



2 

Renaming the variables as u = sin 2 ( 6>3 ~ 6>2 ) and v = sin 2 (y), the above 
equation takes the form 

u 2 (l - v)[(r 2 - R 2 ) 2 + R 2 r 2 uf = v 2 (l - u)[(r 2 - R 2 ) 2 + R 2 r 2 vf. 

This real equation holds only when u — v, that is, sin 2 ( 6 ' 3 ~ 6>2 ) = 
sin 2 (y), or equivalently, 9 3 = 29 2 . If we substitute these values in 
equation of (ITHj) . we obtain 

1 - cos 9 2 (R 4 + r 4 - 2R 2 r 2 cos 29 2 ) 3 



1 - cos 26 2 (R 4 + r 4 - 2R 2 r 2 cos 9 2 f ' 

Taking w = cos #2 and s = cos 29 2) we are led to 

(1 - w)(R 4 + r 4 - 2R 2 r 2 sf = (1 - s)(R 4 + r 4 - 2R 2 r 2 w)\ 

which has real solutions only for w — s, i.e. when cos 9 2 = cos 29 2 , 
which yields 

2ir 47r 

2 = °' T' ¥' 27r ' 

To avoid singular configurations, we must take 9 2 = which 
correspond to an equilateral triangle positioned in the suitable circle of 
radius r. This remark completes the proof. □ 

4. The curved tvbody problem in M 2 r 

Our attempts to study the relative equilibria corresponding to the 
subgroups G\ and G3 in the Poincare disk 3 2 R led to insurmount- 
able computations. Therefore we chose to move the problem to the 
Poincare's upper- half-plane model, M 2 R , to see if the those equilibria 
would be easier to approach. In the remaining part of the paper, we 
will obtain the equations of motion in the Poincare upper half plane 
and analyze the hyperbolic and parabolic relative equilibria for n = 2 
and n = 3. 

4.1. Equations of motion in M 2 R . We will next obtain the equations 
of motion of the curved n-body problem in Poincare's upper half plane 
model, Ef|j, with the help of a global isometric fractional linear trans- 
formation defined by 

/ \ ^9 ^9 , N —Rw + iR 2 

80 z-.Ul^Bl, z = z(w) = — -— , 

w + iR 
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where w is the complex variable in the upper half plane 

m 2 R = {w e c|im(» > o}. 



iR(R-z) 



This transformation has the inverse 

w.Bl^Mi, w = w(z) 

H H w R + z 

Since 

—2R 2 i _ 2R 2 i 

dz = —dw and dz = r^dw, 

(w + 1R) 2 {w - iR) 2 

the metric (Tl3|) of the disk 3 R is transformed into the metric 

(81) — ds 2 = —rdwdw. 

(w — w) z 

Then H|, endowed with the metric (I8T1) is called the Poincare upper 
half plane model of hyperbolic geometry, for which the conformal factor 
is given by 

AR 2 

H(w,w) = — rr^- 

{w — w) z 

In terms of the metric (JHTT) . the geodesies are either half circles orthog- 
onal to the real axis (y = 0) or half lines perpendicular to it. All these 
curves have infinite length. 

Applying transformation (IHUl) to equation (Tl9l) . we obtain the new 
potential in the coordinates (w, w), given by 
(82) 

1 (Wk + w k ){wj + Wj) - 2{\w k \ 2 + \ Wj \ 2 ) 
Vr{w,w) — — > m k mj — — -— , 

where 

(83) Q 3:{kJ )(w,w) 

= {(w k + w k )(wj + wj) - 2(\w k \ 2 + l^jl 2 )] 2 - (w k - w k ) 2 (wj - wj) 2 . 
Notice that with the values w k = x k +iy k and Wj = Xj +iyj, we obtain 

(84) e 3 , (M 

= 4(x k - Xj ) 2 [(x k - Xj) 2 + 2{y 2 k + y 2 )] + 4(y 2 k - y 2 ) 2 , 

which is always positive, except when it takes the value zero as collisions 
take place, i.e. when x k = Xj and y k = yj for at least 2 indices k,j G 
{1, . . . , n}, k 7^ j. 

A result analogous to Lemma [j] in the Poincare disk, 3 2 R , is the 
following property, valid in the Poincare upper half plane, M 2 R . 
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Lemma 3. Let 

1 n 

L(w, w) = 2 mfc ^ Wk ' l^ fc ' 2 + v n( w > ™) 

k=l 

be the Lagrangian for the given problem in Hfj. Then the solution 
curves of the corresponding Euler-Lagrange equations satisfy the system 
of n second order differential equations, k = 1, . . . , n, 

■■ 2m k w\ 2 dV R (w k - w k ) 2 dV R 
85 m k w k — = — — — = — — - , 

w k -w k fi{w k ,w k ) dw k 2R 2 dw k 

where 

fsfil ® Vr - 4m k mjR(w k - w k ){wj - Wj) 2 (w k - Wj)(wj - w k ) 

[ } dw k -jr[ (e 3>{fcj) (^)) 3 / 2 



and G 3) (fej)(w, w) is defined by (JHSj) - 

Proof. From the definition flHUl) of the linear fractional transformation 
z, we obtain in terms of the coordinates that 

-2R 2 i 
AR 2 i 2 2RH 

(w k + iRy (w k + iRy 

dw k _ (w k - iR) 2 
dz k ~ 2R 2 i 

Therefore the equations of motion f[2"51) of the curved n-body problem 
in the Poincare disk get transformed into system (1851) in the Poincare 
upper half plane. This remark completes the proof. □ 

4.2. Integrals of motion in M 2 R . We will further compute the 4 in- 
tegrals of system fl85l) by transforming the 4 integrals of system ()2ip 
via the transformation flHUj) . 
The energy integral is given by 

1 n 

( 87 ) -^2m k fi(w k ,w k ) \ w k \ 2 - V R {w,w) = h, 



k=l 



where V R is given by (1821) and h is the energy constant, the same that 
occurs in (13TI) . which was used to obtain (IB7|) . 
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With the transformation ( )80|) . the 3 integrals of the total angular 
momentum, (1521). (1531. and (EH), become 



Em k R[\w k \ 2 + R 2 -iR{w k - w k )} 2 _ 2 . _ 2 d2/1 . N1 
/on ^ ^m fc i? 2 [|w fc | 2 + # 2 -zi?(w fc -w fc )] 2 / • - , - \ 

(89) g K - + *m - {WkWk + WkWk) = C2 ' 

(9 n 0) 

>A m fc i2[|u;fc| 2 + i2 2 -iR(w k - w k )} 2 „ 2 _ 2 2 ,. „ 



which are the integrals of the total angular momentum for system ( 1851) . 
A straightforward computation confirms that the left hand sides of 
equations f )88|) . fl89|) . and ( 190]) are real functions. Moreover, the con- 
stants ci,C2,C3 G R are the same that occur in equations fl32|) . fl33|) . 
and (El 



4.3. Relative equilibria in Eff^. In this subsection we give conditions 
for the existence of hyperbolic and parabolic relative equilibria in the 
negative curved problem by using the Poincare upper half plane model, 
Elfj. The definitions of these concepts are given in the same geometric 
terms we used in Definitions [T] and [3J 
Let 

SL(2,R) = {Ae GL(2,R) | det A = 1}, 

be the special linear real 2- dimensional group, which is a 3-dimensional 
simply connected, smooth real manifold. It is well known (see, e.g., 
[TT] ) that the group of proper isometries of is the projective quotient 
group SL(2,R)/{±I}. Every class 



A=( a c ^ eSL(2,R)/{±J} 



has also associated a unique Mobius transformation ] a '■ HPjj — > M 2 R , 
where 

az + b 

M*) = — Zl> 

cz + a 

for which it is easy to see that /-a( z ) = Ia( z )- The Lie algebra of 
SL(2,R) is the 3-dimensional real linear space 

sl(2, R) = {X G M(2, R) | trace X = 0}, 
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spanned by the following suitable set of Killing vector fields, 

*1 = \ ( l _°1 J > X2 = ( j ) ' X 3 = ( -1 l 

As in section 13. 3^ we consider the exponential map of matrices, 

exp : *Z(2,R) -> SL(2,R), 

applied to the one-parameter additive subgroups (straight lines) {tXi}, 
{tX 2 }, and {tX 3 }, to obtain the following one-parameter subgroups of 
the Lie group SL(2,R): 

(1) The isometric dilatation subgroup 

/ e t/2 
0i (t) = exp(tXi) = ( o e -t/2 

which defines the one-parameter family of acting Mobius trans- 
formations 

(91) f 1 (w,t)=e t w(t); 

(2) The isometric shift subgroup 



(h{t) = exp{tX 2 ) 



1 t 
1 



which defines the one-parameter family of acting Mobius trans- 
formations 

(92) f 2 (w,t) = w(t) + t; 

(3) The isometric rotation subgroup 



03 (t) = exp(tX 3 



cos t sin t 
- sin t cos t 



which defines the one-parameter family of acting Mobius trans- 
formations 

(cost) w(t) + sin t 

(93) f 3 (w,t) = 



— sint) w(t) + cost 

We proved in subsection 13 .41 the existence of elliptic relative equilibria 
for the initial problem by using the proprieties of the Poincare model 
3 2 R , i.e. we showed that for elliptic relative equilibria, each particle 
moves along a circle centered at the origin of the coordinate system. 

From the Theorem of the invariance of the domain, [H], the isom- 
etry ( |80|) carries the interior of into the interior of M 2 R . Therefore 
simple closed curves contained in B) R are taken to simple closed curves 
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contained in IHrJj (see [2] for more details). For the circle z(t) = zoe lt 
in the Poincare disk, 3 R , the corresponding curve, 

, . iR(R - z e lt ) 

w(t) = ■ 

in the Poincare upper half plane, M R , satisfies w(0) = w(2tt) and must 
therefore belong to the class of the Mobius transformations (|93|) . which 
thus corresponds to elliptic relative equilibria. Since we already studied 
those orbits in the Poincare disk, D^, we don't need to further analyze 
them here, the more so since the Mobius transformations (1931) lead to 
complicated computations in M R . However, the Mobius transforma- 
tions (191 p and (|9~2j) . corresponding to the Killing vector fields Xi and 
X2, respectively, are simpler in M 2 R than the analogue transformations 
in 3 2 R , so they will be the object of our further analysis. The former 
will lead us to hyperbolic relative equilibria and the latter to parabolic 
relative equilibria. 

4.4. Hyperbolic relative equilibria. We will further study the rel- 
ative equilibria associated to the subgroup (}9Tj) . which defines the one- 
parameter family of acting Mobius transformations 

f 1 (w,t) = e t w(t) 

in the Poincare upper half plane, H^. Let £ = (Ci, ■ ■ ■ ,Cn), with 
£k(t) = e t w k {t), k = 1, . . . ,n, be the action orbit for a solution w = 
(wi, . . . , w n ) of system ( 1851) . Then, 

£fc = (w + w)e\ £ fc = (w + 2w + w)e\ k — 1, . . . , n, 

and therefore the curve £ is also a solution of the equations of motion 
(|85|) if and only if 

2m k e 2t (w k + w k ) 2 (e l w k - e l w k ) 2 dV R dw k 



m k (w k + 2w k + w k )e 

2m k (w k + 2w k w k + wl)e l (w k - w k ) 2 dV R t 



e l w k - e l w k 2R 2 dw k d£ k 

k = 1, . . . , n. 



w k — w k 2R 2 dw k 

Since w is a solution of fl85l) . — =^ = e~ l 7^ 0, and m k ^ 0, k — 1, . . . , n, 
the last relationships become 

(94) 2w fc — = = w k , k = l,...,n. 

w k ~w k w k - w k 

So if we fix a body m k , the above condition holds if and only if 

2w k 

1 — = 

w k - w k 
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or 

(95) 2w k = -w k . 

The first condition in equation is equivalent to w k + w k = and 
corresponds to the geodesic given by the vertical half line that forms 
the imaginary axis. Via the linear fractional transformation flHUl) . this 
geodesic corresponds to the horizontal geodesic, Im(2;) = 0, of the 
Poincare disk B) R . 

The second condition, fl95|) . holds for the body m k when w k is a 
function of the form w k (t) = W7£,.(0)e~*/ 2 , where Wk(0) is some initial 
condition with lm(wk{0)) > 0. Therefore the particle m k moves along a 
half line through the origin of the coordinate system, which is a point 
at infinity that is reached when t — > oo. As t — > — oo, m k goes to 
infinity. It is instructive to notice that such integral curves of equation 
fl53|) correspond in the Poincare disk to the parametric curves 

-Rw k (0)e^ 2 + iR 2 _ iR 2 e l l 2 - Rw k (0) 
Zk ^> ~ ^(0)e-*/ 2 + iR ~ iRe^ + WkiO) ' 

which are equidistant from the horizontal geodesic diameter, all of them 
starting at the point (— R, 0), as t — > — oo, and ending in the point 
(R, 0), as t — > oo. When Re[wfc(0)] = 0, this curve becomes the hor- 
izontal diameter, which is a geodesic. The above non-geodesic curves 
together with the geodesic horizontal diameter foliate the Poincare disk. 

On the hyperbolic sphere, h 2 R , the above non-geodesic curves are 
given by the equations 

(96) r = a, rj = rsinht, 3 = rcosht, 



where a ^ is a constant and r = yj R 2 + a 2 . When a = 0, we recover 
the geodesic. When these curves are taken to the Poincare disk via the 
stereographic projection (j3J), we obtain the complex curves 

2Rr(e 2t -l) Rae 1 i 

(97) At) = n { t) + ,„ W = - - I T) + - - 2r(e2 , - iy 

which also start at the point (— R, 0), as t — > — oo, and end at the 
point (R, 0), as t — > oo. In other words, we have obtained two families 
of equidistant curves with the same initial and final directions as the 
geodesic diameter. 

So for the particles mi, . . . , m n to form a relative equilibrium in M 2 R 
associated to the Killing vector field X%, they have to move along the 
upper half lines converging to the origin of the coordinate system in 
Hfjj as t — > oo. These half lines are not geodesies, except in the case of 
the vertical half line. Moreover, each non-geodesic vertical half line is 
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equidistant from the geodesic vertical half line, the distance being larger 
when the angle between the non-geodesic half line and the geodesic 
vertical half line is larger. The sizes of these angles range between 
and 7r/2. Notice also that in the case of the geodesic vertical half line, 
equation ( 195]) is satisfied as well, a fact that we will use in the proof of 
the following result. 

Theorem 6. Consider n point particles with masses mi, ... , m n > 0, 
n > 2, moving in M 2 R . A necessary and sufficient condition for the 
function w = (wi, . . . ,w n ) to be a solution of system (1H5|) and, at 
the same time, a relative equilibrium associated to the Killing vector 
field Xi defined by equation (19 ip is that, for every k = 1, . . . ,n, the 
coordinates satisfy the conditions 

R(w k + w k ) w k _ ^ rnjjwj - Wj) 2 (w k - Wj)(wj - w k ) 

where 

(99) e 3i{kJ) (w,w) 

= [(w k + w k )(wj + wj) - 2(\w k \ 2 + |^j| 2 )] 2 - (w k - w k ) 2 (wj - Wj) 2 , 
k,j e {1, • • • ,n}, j. 

Proof. We showed previously that for relative equilibria of the afore- 
mentioned type the bodies move along straight half lines converging to 
the origin of the coordinate system and must therefore satisfy equation 
(|95|) . which implies that 

(100) 4w k = w k , k = l,...,n. 

Using this equation together with equation (195|) . we can conclude from 
the equations of motion (|85l) that 

nm x R 2 m k (w k + w k ) w k _ dV R 

{WL) — — — - — — — , K— L,...,n. 

2(w k - Wk) 6 dw k 

Using ( |86i) . we obtain the relationships given in ( 1981) . This remark 
completes the proof. □ 

Definition 3. We will call hyperbolic relative equilibria the solutions 
of system (!85l) in H|, that satisfy equations (T98T) . 

We remark that equation (1951) also gives the condition the velocity of 
the particle m k must satisfy in order to produce a hyperbolic relative 
equilibrium. 
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4.4.1. The case n = 2. We will further provide a description of the 
hyperbolic relative equilibria for 2 interacting particles in M 2 R . For this, 
we observe that for the particles of masses mi and m 2 , the equations 
which characterize hyperbolic relative equilibria become 

R(wi + wx) wi m 2 (w 2 - w 2 ) 2 (wi - w 2 )(w 2 - w x ) 



(102) 



(103) 



8K- Wl ) 4 [e 3 ,(i, 2) (^)] 3 / 2 

R(w 2 + w 2 ) w 2 mi(u>i - wi) 2 (w 2 - wi)(m)i - w 2 ) 



8(w 2 -w 2 y [e 3 ,( 2 ,i)(^)] 3 / 2 
where 

®3,(i,2)(w,w) = @3,(2,i)(w,u;) = A(w 2 -w 1 )(wi-w 2 )(wi-w 2 )(w 2 -wi). 
Straightforward computations lead to the equation 

^ 1Q4 ^ (wi + wi)wi _ m 2 (wi - wi) 2 (w 2 - wi) 

(w 2 + w 2 )w 2 mi(w 2 - w 2 ) 2 (wi - w 2 ) ' 
provided that w 2 + w 2 ^ 0, which is equivalent with the equation 
(105) m 1 (w 2 - w 2 ) 2 (wi - w 2 )(wi + wi)w\ 

+m 2 (wi - wi) 2 (w 2 - wi)(w 2 + w 2 )w 2 = 0, 
subject to the restriction w 2 + w 2 ^ 0. 

Let us first prove a negative result about hyperbolic relative equilib- 
ria in the curved 2-body problem in Hf|j. In general terms, unrelated 
to any model of hyperbolic geometry, this result states that, on one 
hand, 2 particles cannot follow each other along a geodesic and main- 
tain a constant distance between each other and, on the other hand, 
one particle cannot move along the geodesic, while the other particle 
moves along a non-geodesic curve equidistant from that geodesic, such 
that the bodies maintain all the time the same distance between each 
other. 

Proposition 3. Consider 2 point particles of masses m 1; m 2 > mov- 
ing in EI^. Then there are no hyperbolic relative equilibria as solutions 
of system ( 1851) with n = 2 for which both particles move along the geo- 
desic vertical half line on the imaginary axis or for which one particle 
moves along the geodesic vertical half line and the other particle moves 
along a non-geodesic half line converging to the origin of the coordinate 
system. 

Proof. As previously seen, if the particle of mass mi moves along ge- 
odesic vertical half line, then wi + Wi = 0. Therefore equation f )105p 
becomes 

m 2 (wi - wi) 2 (w 2 - wi)(u> 2 + w 2 ) w 2 = 0, 
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which is impossible because none of the above factors vanishes in either 
of the two scenarios given in the statement. This remark completes the 
proof. □ 

The previous result showed that 2 bodies that form a hyperbolic 
relative equilibrium cannot move along the same geodesic. It is then 
natural to ask whether they could move along the same non-geodesic 
curve that is equidistant from a given geodesic and maintain all the 
time the same distance from each other. As we show in the following 
result, expressed in terms of M 2 R , the answer is also negative. 

Proposition 4. Consider 2 point particles of masses m 1; m 2 > mov- 
ing in Hjj. Then there are no hyperbolic relative equilibria as solutions 
of system (185j) with n = 2 for which both particles move along the same 
half line converging to the origin of the coordinate system. 

Proof. Recall that for hyperbolic relative equilibria the expression of 
the orbit of the particle is Wk(t) = u>fc(0)e~'/ 2 , where lm(wfc(0)) > 0. 
Then, if 2 particles belong to the same half line and form a hyperbolic 
relative equilibrium, their corresponding solutions have the form 

(106) loi(f) = wi(0)e~* /2 and w 2 (t) = aw x (t) = awi(0)e" t/2 , 

for some a > 0. If we substitute relations (I106P into equation (11051) . a 
straightforward computation gives us the condition 

(ami + m 2 )u>i(0) = (mi + am 2 )tui(0). 

If Re(wi(0)) 7^ 0, i.e. when the particles are on the same non-geodesic 
half line, then there is no a that can satisfy the equation. This remark 
completes the proof. If Re(wi(0)) = 0, then the bodies are on the 
geodesic vertical half line, and only a = — 1 satisfies the equation. But 
since a must be positive, the equation is not satisfied either, so we 
found another proof for the first statement of Proposition [3l □ 

From Propositions [3] and H] we learned that, in general terms, inde- 
pendent of the hyperbolic model used, 2 particles cannot form a hy- 
perbolic relative equilibrium if they move along the same non-geodesic 
curve equidistantly places from a geodesic, along the same geodesic, or 
if one particle moves on a geodesic and the other particle on a non- 
geodesic curve equidistant to the geodesic. It is then natural to check 
the last possibility, whether there exist hyperbolic relative equilibria 
with one body moving along a non-geodesic curve and the other along 
another non-geodesic curve, both equidistant from a given geodesic. 
The answer is positive and the motion takes place when the distances 
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from the geodesic to the 2 non-geodesic curves satisfy a certain rela- 
tionship that depends on the values of the masses. 

To write these conditions in the language of M 2 R , recall that the slope 
/3 of any straight line in the complex plane is defined by the formula 

n w — w 

ip = -. 

w + w 

Without loss of generality, we choose the initial conditions such that 
the heights satisfy wi(0) — u>i(0) = 2iyi and u>2(0) — ^2(0) = 2iy 2 . 
Then, in terms of the slopes /3i and (3 2 of the straight lines, equation 
(I105P becomes 

(107) mi f3 2 y 2 {\wi\ 2 - w 2 w 1 ) = -m 2 /3i yi (\w 2 \ 2 - w 2 wi). 
We can now state and prove the following result. 

Theorem 7. Consider 2 point particles of masses m 1; m 2 > moving 
in H|j. Then some necessary and sufficient conditions for the existence 
of a hyperbolic relative equilibrium as a solution of system ( 1851) with 
n = 2 are that one particle moves along a non-geodesic half line, while 
the other particle moves along another non-geodesic half line, both half 
lines converging to the origin of the coordinates system, such that the 
supporting lines have slopes of opposite signs that satisfy the relation- 
ship 

(108) ^ = -^ y A 

m 2 (3i y 1 

and that at every time instant there is a geodesic half circle centered at 
the origin of the coordinate system on which both particles are located. 

Proof. If we take the real and imaginary parts of equation f)107p . we 

obtain the equations 

(109) 

"ii/%2|u>i| 2 - m 1 (3 2 y 2 Re(w 2 w 1 ) = -m^y^w^ 2 + m 2 (3 1 y 1 Re(w 2 w 1 ), 

(110) - mi/3 2 2/2lm(w2^i) = , m 2 /3iy 1 Im(w 2 wi). 

From equation fill 01) . it follows that 

-m 1 (3 2 y 2 = m 2 f3iy 1 , 
so condition (11081) must be satisfied. Equation fll(J9j) implies that 

rai/%2M 2 = -m-i^y^w^ 2 , 

which, by condition (11081) . is equivalent to |wi| 2 = \w 2 \ 2 . This fact 
proves that, at every time instant, there must exist a geodesic half 
circle centered at the origin of the coordinate system on which the 
particles are located. This remark completes the proof. □ 
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Remark 2. Notice that, in the above result, if the masses are equal, 
m i — m 2 > 0, then both half lines have slopes equal in absolute value, 
but of opposite sign. Indeed, since, at every time instant, the particles 
must be located on the same geodesic half circle centered at the origin 
of the coordinate system, the slopes fa and fa of the half lines along 
which they move satisfy Pi = — fa if and only if y\ = y 2 . Therefore, 
from equation (I108p . we have that ni\ = m 2 if and only if fa = —fa. 

4.4.2. The case n = 3. We will next study the case of 3 bodies in the 
Poincare upper half plane, H^, with masses mi,m 2 ,m 3 > 0. In this 
context, the system of algebraic equations becomes 

(111) R^ + wJw, 

m 2 (w 2 - w 2 ) 2 {wi - w 2 )[w 2 - wi) 
{[(wi + wi)(w 2 + w 2 ) - 2(\w 2 \ 2 + \wi\ 2 )} 2 - (u>i - wi) 2 (w 2 - w 2 ) 2 } 3 / 2 

m 3 (w 3 - w 3 ) 2 (wi - w 3 )(w 3 - wi) 

{[(wi + wi)(w 3 + w 3 ) - 2(\w 3 \ 2 + |wi| 2 )] 2 - (wx - wi) 2 (w 3 - ms) 2 } 3 / 2 ' 

R(w 2 + w 2 ) w 2 



(112) 



8(w 2 - w 2 y 

m\{wi - wi) 2 [w 2 - wx){wi - w 2 ) 



{[{w 1 + wi)(w 2 + w 2 ) - 2(|u; 2 | 2 + l^il 2 )] 2 - (wi - wi) 2 (w 2 - u; 2 ) 2 } 3/2 

m 3 (w 3 - w 3 ) 2 {w 2 - w 3 )(w 3 - w 2 ) 

{[(w 2 + w 2 )(w 3 + w 3 ) - 2(\w 2 \ 2 + \w 3 \ 2 )] 2 - (w 2 - w 2 ) 2 (w 3 - w 3 ) 2 } 3 / 2 ' 

R(w 3 + w 3 ) w 3 



(113) 



8(w 3 - w 3 ) 4 
m x {wi - wi) 2 {w 3 - w\){wi - w 3 ) 



{[(wi + wi)(w 3 + w 3 ) - 2(\w 1 \ 2 + |^ 3 | 2 )] 2 - (wi - wi) 2 (w 3 - u; 3 ) 2 } 3 / 2 
m 2 (w 2 - w 2 ) 2 (w 3 - w 2 )(w 2 - w 3 ) 



{[(w 2 + w 2 )(w 3 + w 3 ) - 2(\w 2 \ 2 + |w 3 | 2 )] 2 - (w 2 - w 2 ) 2 (w 3 - w 3 ) 2 } 3 / 2 
We assume that the particles mi and m 3 move along non-geodesic 
half lines and that m 2 moves along the geodesic vertical half line, such 
that, at every time instant, there is a geodesic half circle centered at 
the orgin of the coordinate system on which all particles are located. 
In other words, if w = (wi,w 2 ,w 3 ) represents the configuration of the 
system, we have 

w 2 + w 2 = and \wi\ = lii^l = \ w 3\- 
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We can also write that 

Mt) = ^i(0)e-* /2 , w 2 (t) = w 2 (0)e- t/2 , w 3 (t) = w 3 (0)e~ t/2 , 

with |wi(0)| = |u> 3 (0)| = 1 and w 2 (0) = i- The latter conditions are not 
restrictive, since the only requirement for the initial conditions is to lie 
on the half lines on which the particles are assumed to move. Substi- 
tuting the above forms of wi,w 2 ,w 3 into equations ( 11111) . ( 11121) . ( 11131) . 
the factors e - */ 2 get cancelled, and after redenoting u>i(0), ^2(0), u» 3 (0) 
by iwi, w 2 , w 3 , respectively, we obtain the following equations: 

R(wi + wi) wi m 2 (w 2 + l) 



(114) 



;(wi-wi) 4 ~ 2[A + (w 1 -w 1 ) 2 ] 3 / 2 

m 3 {w 3 - w 3 ) 2 Oi - w 3 )(w 3 - w x ) 



(115) 
(116) 



{[(w x + w 1 )(w 3 + w 3 ) - 4] 2 - (wi - Wi) 2 (w 3 - w?) 2 }*/ 2 ' 

_ mi(wi - wi) 2 (wi + wi) m 3 (w 3 - w 3 ) 2 (w 3 + w 3 ) 
[4 + ( Wl - w^ 2 ] 3 / 2 [4 + (w 3 - w 3 ) 2 Y' 2 

R(w 3 + w 3 ) w 3 _ m 2 (wl + 1) 
8(w 3 -w 3 )* ~ 2[A+{w 3 -w 3 ) 2 fl 2 
m 1 (w 1 -wi) 2 (w 3 -wi)(w 1 -w 3 ) 



{[{wi + w 1 )(w 3 + w 3 ) - 4] 2 - (wi - wi) 2 (w 3 - w 3 ) 2 } 3 / 2 ' 
Equation (11151) can take place only if Wi + Wi and w 3 + w 3 have 
opposite signs, so we can rewrite this equation as 

m 1 (w 1 - wi) 2 |wi + wi\ _ m 3 (w 3 - w 3 ) 2 |w 3 + w 3 \ 

[ } [2K+»i)F ~ [2(w 3 + w 3 )] 3 

To express w\ and w 3 in terms of the angles the half lines make with 
the horizontal axis, we put 

Wi = e ldl = cos Q\ + % sin 6\ and w 2 = e ld2 = cos 9 2 + i sin 9 2 , 

with 9 U 9 3 e (-tt/2,0) U (0,tt/2). Then equation (TTTT|) becomes 

(118) mi tan 2 9i = m 3 tan 2 9 3 , 

which shows what relationship exists between the masses and the angles 
of the non-geodesic half lines along which the corresponding particles 
move. 

We can now state and prove the main result of this section, which 
shows that Eulerian relative equilibria for which one body moves along 
a geodesic exist only if the masses moving on non-geodesic curves 
equidistant from the geodesic are equal and those curves are on op- 
posite parts of the geodesic and at the same distance from it. 
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Theorem 8. Consider 3 point particles of masses 7711,7712,7713 > 
moving in H^. Assume that mi and m 3 move along non-geodesic half 
lines emerging from the origin of the coordinate system at angles Q\ and 
83, respectively, and thatm 2 moves along the geodesic vertical half line. 
Moreover, at every time instant, there is a geodesic half circle on which 
all 3 bodies are located, and the motion of the particles is given by the 
function w = (wi, 102,103). Then w is a hyperbolic relative equilibrium 
that is a solution of system (185jl with n = 3 if and only if Q\ = —63 
and mi = m 3) with 9 U 9 3 £ (-7r/2, 0) U (0, 7r/2). 

Proof. We already showed that for wi = e l9l ,w 2 = i, and w 3 = e^ 3 , 
equation (11151) takes the form fll 18j) . With the same substitutions, 
equations fll 14|) and fll 16j) become, respectively, 

. . Rcos9i 8m 2 Smssin 2 ^ 

( 119 ) _ : _4/i = „„„2/i + 



sin 9i cos 2 61 (cos ^3 — cos6 l i) : 
R cos #3 8mi sin 2 9i 8m 2 



^ ^ sin 4 6*3 (cos 9i — cos 6*3)2 ' cos 2 #3 

If we divide equation fll 19j) by cos 2 #3 and equation f ll 2Q[) by cos 2 #i, 
using relation fl 1 1 8 j) we obtain that 

COS 3 gi _ COS 3 fl 3 
v / -4/1 ~~ ■ 4/i ■ 

sin t^i sm #3 

Consider now the function 



/: (-vr/2,0)U(0,vr/2), f{x) 



COS 3 X 
~. 4 : 

sm x 



which is obviously even. It is easy to see that / is increasing in the 
interval (— 7r/2,0) and decreasing in the interval (0, 7r/2). Therefore 
equation f 1 1 2 1 j) has solutions if and only if 9i = ±#3. Since #1 = #3 
induces a collision configuration, which is a singularity, the only pos- 
sible solution is 9\ = —#3. The fact that mi = m 3 follows now from 
equation (11181) . This remark completes the proof. □ 

4.5. Parabolic relative equilibria. In this section we will study the 
relative equilibria associated to the subgroup 



fc(t) = ex V (tX 2 ) 



1 t 
1 



generated by the Killing vector field X 2 and which defines the one- 
parametric family of acting Mobius transformations 

f 2 {w,t)=w{t)+t, 



40 



F. Diacu, E. Perez-Chavela, and J.G. Reyes Victoria 



in the upper half plane M 2 R . These orbits correspond to parabolic rel- 
ative equilibria, and we will show that they do not exist in H|,. 

Let C — (Ci) • • ■ > Cn)> with Ck(t) = Wkif) + t, be the action orbit for 
w = (wx, . . . , w n ), which is a solution of the equations of motion (!85|) . 
Then 

Cfe = ™k + 1 and (k = w k , k = 1, . . . , n, 
therefore ( is also a solution of system (1851) if and only if 

mtG = = n — , k = 1, • • • , n, 

which can be written as 

2m k (w k + l) 2 - w h ) 2 dV R dw k 
rrikWk = z Tvi 77~= ^> k = 1, . . . , n, 

which, since = 1, is the same as 
d(k 

n nnN .. 2m fc (Wfc + 2w fc + 1) O fe - w fc ) 2 9Vr 

(122) m fc w; fc = — — - , k = 1, . . . , n. 

w k - w k 2R A dw k 

But since w is a solution of system fl85l) , we also have that 

(123) mA= i!^-(«)!^, fc= !,...,„, 

Comparing now equations (11221) and (11231) . we obtain that 

(124) 2w k = -1, fc = 1, . . .,n, 
which holds if and only if 

(125) w k {t) = -- + w fc (0), fc = 1, . . . , n, 

where Wk{0), k = 1, . . . , n, are initial conditions. Consequently, a nec- 
essary condition for the particles m 1; . . . , m n to form a relative equilib- 
rium associated to the Killing vector field X 2 is that they move along 
horizontal straight lines in M 2 R passing through Wk(0), k — 1, . . . , n. In 
terms of the Poincare disk B) R , these orbits correspond to the paramet- 
ric curves 

Zk[t) - -* +Wk (Q) + iR ^-^■■■^ 

which start at the point (— R, 0), as t — > — oo, and end at the same 
point (-R, 0), as t — > oo. These curves have the same topology as the 
boundary circle of 3 R . In terms of the hyperbolic sphere h 2 R , these lines 
are the parabolas obtained by intersecting L|, with a plane orthogonal 
to the rotation axis, the line y = 0, z — x. 
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We can now state and prove the following result. 

Theorem 9. Consider n > 2 point particles of masses m 1; . . . , m n > 
moving in M 2 R . Then a necessary and sufficient condition for the 
function w = (wi, . . . ,w n ) to be a solution of system fl85|) that is a 
relative equilibrium associated to the Killing vector field X 2 is that the 
coordinate functions satisfy the equations 

/ 12g N R _ m j{wj ~ w j ) 2 {w k - Wj)(wj - w k ) 

~ ±(w k -w k y~ [e 3i(M ( w ,^)] 3 / 2 

k — 1, . . . , n, where 

(127) 03,(fe,i) (w,w) = [{w k + w k )(wj + Wj) - 2{\w k \ 2 + |^j| 2 )] 2 

-(w k - w k ) 2 (wj - wj) 2 , k,j E {1, . . . ,n}, k^j. 

Proof. We saw that relative equilibria associated with the Killing vector 
field X 2 must satify equations (I125p . which imply that 

(128) w k = 0, k — l,...,n. 

Therefore, from equations (11231) . we can conclude that the coordinates 
of a relative equilibrium w satisfy the equations 

m k R 2 8Vr 
129 k —- = —^, fc = l,...,n. 

[w k - w k ) 6 dw k 

If we now compare the above equations to the expressions (1861) of 
l^p, k = l,...,n, we obtain the desired relationships (11271) . This 
remark completes the proof. □ 

Definition 4. We will call parabolic relative equilibria the solutions of 
system (I85p in M 2 R that satisfy equations 



Notice that equations (11241) provide the velocities of the particles in 
case they form a parabolic relative equilibrium. However, as we will 
further prove, parabolic relative equilibria do not exist in the curved 
77,-body problem. The following statement generalizes a result obtained 
in [8] for curvature k — — 1. Using the same idea as in [8], this result 
was generalized in [3] to the 3-dimensional case. 

Theorem 10. In the curved n-body problem with negative curvature 
there are no parabolic relative equilibria. 

Proof. Using the notation 

w k = a k + ib k and Wj = Oj + ibj, 
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the real and imaginary part of equations (11261) become 
(m) R _^ 2m ] b%a k -a ] f + {b)-bl)] 

326 4 fc [e 3 , {kj) (w ,w )]^ 



n, 



^ 131 j Q =J2 Amjb2 ^ ak ~ aj 



J^l [@3,(fcj)(w ,W )] 3/2 



n. 



Since bk > 0, k = l,...,n, equations fl 1 3 1 [) hold for any k,j G 
{1, . . . , n}, with k ^ j, if and only if a k = aj. This fact implies that 
all the particles are located on the same vertical line. Without loss 
of generality, we can assume that they are on the vertical half line 
x = 0, y > 0. Therefore w k = b k i and Wj = bji. When we substitute 
these values into equation (11261) we obtain 

(132) R -f ^glj k-1 n 

{ } 32^"^ \bl-b^ ' k - l --- n - 

Since the particles do not collide, we can assume, without loss of gen- 
erality, that 

< &i < • • • < b n . 
Then, for k = 1, we can conclude from (11321) that 

R _™ m 3 b){b)-b\) 
32bi~^ \b\-%\* ■ 

But the left hand side of this equation is negative, whereas the right 
hand side is positive. This contradiction completes the proof. □ 
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